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Abstract. Hubble Space Telescope (HST)[] colour - mag- 
nitude diagrams in B, V and R along with long-slit Mul- 
tiple Mirror Telescope (MMT)Q spectrophotometric data 
are used to investigate the evolutionary status of the 
nearby blue compact dwarf (BCD) galaxy I Zw 18. We 
find that the distance to I Zw 18 should be as high as 20 
Mpc, twice the previously accepted distance, to be con- 
sistent with existing observational data on the galaxy: 
colour-magnitude diagrams, the high ionization state of 
the gas and presence of WR stars in the main body, and 
the ionization state of the C component. The spectral en- 
ergy distribution (SED) of the main body of I Zw 18 is 
consistent with that of a stellar population with age Si 5 
Myr. However, the presence of large-scale shells observed 
around the main body suggests that star formation began 
~ 20 Myr at the NW end and propagated in the SE di- 
rection. Our analysis of colour-magnitude diagrams and 
of the spectral energy distribution of the C component 
implies that star formation in this component started ^ 
100 Myr ago at the NW end, propagated to the SE and 
stopped ~ 15 Myr ago. Thus, I Zw 18 is likely to be one 
of the youngest nearby extragalactic objects. 

Key words: galaxies: evolution - galaxies: formation - 
galaxies: ISM - galaxies: photometry - galaxies: starburst 
- galaxies: stellar content 



1. Introduction 

I Zw 18 remains the most metal-poor blue compact dwarf 
(BCD) galaxy known since its discovery by Sargent & 
Searle (1970). Later spectroscopic observations by Searle 

Send offprint requests to: Y. Izotov 

1 Based on observations obtained with the NASA/ESA Hub- 
ble Space Telescope through the Space Telescope Science Insti- 
tute, which is operated by AURA, Inc. under NASA contract 
NAS5-26555. 

2 Ground-based spectroscopic observations presented herein 
were obtained with the Multiple Mirror Telescope, a facility 
operated jointly by the Smithsonian Institution and the Uni- 
versity of Arizona. 



& Sargent (1972), Lequeux et al. (1979), French (1980), 
Kinman & Davidson (1981), Pagel et al. (1992), Skill- 
man & Kennicutt (1993), Martin (1996), Izotov, Thuan 
& Lipovetsky (1997c), Izotov & Thuan (1998), Vi'lchez & 
Iglesias-Paramo (1998), Izotov & Thuan (1999) and Izo- 
tov et al. (1999) have confirmed its low metallicity with 
an oxygen abundance of only ~ 1/50 the solar value. 

Zwicky (1966) described I Zw 18 as a double system 
of compact galaxies, which are in fact two bright knots of 
star formation with an angular separation of 5'.'8. These 
two star-forming regions are referred to as the brighter 
northwest (NW) and fainter southeast (SE) components 
(Fig. |l|) and form what we will refer to as the main body. 
Later studies by Davidson, Kinman & Friedman (1989) 
and Dufour & Hester (1990) have revealed a more com- 
plex optical morphology. The most prominent diffuse fea- 
ture, hereafter component C (Fig. Q), is a blue irregular 
star-forming region ~ 22" northwest of the NW compo- 
nent. Dufour, Esteban & Castaneda (1996a), van Zee et 
al. (1998) and Izotov & Thuan (1998) have shown the C 
component to have a systemic radial velocity equal to that 
of the ionized gas in the NW and SE components, thus es- 
tablishing its physical association to I Zw 18. Furthermore, 
van Zee et al. (1998) have shown that this component is 
embedded in a common H I envelope with the main body. 

Searle & Sargent (1972) and Hunter & Thronson 
(1995) have suggested that I Zw 18 may be a young 
galaxy, recently undergoing its first burst of star forma- 
tion. The latter authors concluded from HST images that 
the colours of the diffuse unresolved component surround- 
ing the SE and NW regions are consistent with a popula- 
tion of B and early A stars, i.e. with no evidence for older 
stars. 

Ongoing star formation in the main body of I Zw 18 
is implied by the discovery of a population of Wolf-Rayet 
stars in the NW component (Izotov et al. 1997a; Legrand 
et al. 1997; de Mello et al. 1998). Flux-calibrated opti- 
cal spectra of the C component (Izotov & Thuan 1998; 
van Zee et al. 1998) reveal a blue continuum with weak 
Balmer absorption features and faint Ha and H/3 in emis- 
sion. Such spectral features imply that the HII region is 
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ionized by a population of early B stars and suggest that 
the C component is older than the main body of I Zw 18. 
Izotov & Thuan (1998) have suggested an age sequence 
from the C component ( ~ 200 Myr ) to the SE region ( 
~ 5 Myr ) of active star formation. 

Dufour et al. (1996b) have discussed new HST imagery 
of I Zw 18, including the C component, which is resolved 
into stars. Based on the analyses of colour-magnitude di- 
agrams, they concluded that star formation in the main 
body began at least 30 - 50 Myr ago and is maintained to 
the present, as is apparent in the SE component. Martin 
(1996) and Dufour et al. (1996b) have discussed the prop- 
erties of expanding superbubbles of ionized gas driven by 
supernova explosions and have inferred dynamical ages of 
respectively 15-27 Myr and 13-15 Myr. As for the 
age of the C component, Dufour et al. (1996b) found in 
a (B — V) vs. V colour-magnitude analysis a well defined 
upper stellar main sequence indicating an age of the blue 
stars of ~ 40 Myr. However, numerous faint red stars were 
also present in the colour-magnitude diagram implying an 
age of 100 - 300 Myr. Dufour et al. (1996b) therefore con- 
cluded that the C component consists of an older stellar 
population with an age of several hundred Myr, but which 
has experienced recently a modest starburst in its south- 
eastern half as evidenced by the presence of blue stars and 
Ha emission. 

Recently Aloisi, Tosi & Greggio (1999) have discussed 
the star formation history in I Zw 18 using the same HST 
WFPC2 archival data ( i.e. those by Hunter & Thron- 
son (1995) and Dufour et al. (1996b)). They compared 
observed colour-magnitude diagrams and luminosity func- 
tions with synthetic ones and concluded that there were 
two episodes of star formation in the main body, a first 
episode occuring over the last 0.5 - 1 Gyr, an age more 
than 10 times larger than that derived by Dufour et al. 
(1996b), and a second episode with more intense activity 
taking place between 15 and 20 Myr ago. No star forma- 
tion has occurred within the last 15 Myr. For the C com- 
ponent, Aloisi et al. (1999) estimated an age not exceeding 
0.2 Gyr. 

Garnett et al. (1997) have presented measurements of 
the gas-phase C/O abundance ratio in both NW and SE 
components, based on ultraviolet spectroscopy with the 
Faint Object Spectrograph (FOS) onboard HST. They de- 
termined values of log C/O = -0.63 ± 0.10 for the NW 
component and log C/O = -0.56 ± 0.09 for the SE compo- 
nent. These ratios, being significantly higher than in other 
metal-poor irregular galaxies, apparently require that car- 
bon in I Zw 18 has been enriched by an older generation 
of stars. Garnett et al. (1997) concluded that I Zw 18 
must have undergone an episode of star formation several 
hundred million years ago. 

Using the same HST ultraviolet spectra, Izotov & 
Thuan (1999) have rederived the C/O abundance ratio 
in both NW and SE components of I Zw 18. They ob- 
tained lower values of log C/O equal to -0.77 ± 0.10 for 



the NW component and -0.74 ± 0.09 for the SE com- 
ponent. With these lower C/O ratios, I Zw 18 does not 
stand apart anymore from other low-mctallicity BCDs. 
Furthermore, the C/O ratios are in excellent agreement 
with those predicted by massive star nucleosynthesis the- 
ory. Therefore, no preceding low-mass carbon-producing 
stellar population needs to be invoked, thus, supporting 
the original idea that I Zw 18 is a young galaxy undergo- 
ing its first episode of star formation (Searle & Sargent 
1972). The main source of the differences (~ 0.2 dcx) 
with the values derived by Garnett et al. (1997) comes 
from the adopted electron temperatures. Izotov & Thuan 
(1999) use higher electron temperatures (by 1900 K and 
2300 K respectively for the NW and SE components) as 
derived from recent MMT spectral observations in aper- 
tures which match more closely those of the HST FOS 
observations used to obtain carbon abundances. 

Thus, despite the extensive multiwavelength studies 
of I Zw 18, its evolutionary status remains controversial. 
Increasing evidence is accumulating, however, in favor of 
the idea that this BCD underwent its first episode of star 
formation less than 100 Myr ago. The first line of evi- 
dence is based on heavy element abundance ratios. Izo- 
tov & Thuan (1999) have studied these ratios in a sam- 
ple of low-metallicity BCDs. They found that all galaxies 
with heavy element mass fraction Z Ja Z Q /20, including 
I Zw 18, show constant C/O and N/O abundance ratios 
which can be explained by element production in massive 
stars ( M ^ 9 M Q ) only. Intermediate-mass stars ( 3 
Mq Si M < 9 M Q ) in these galaxies have not had time 
to die and release their C and N production. Izotov & 
Thuan (1999) put an upper limit to the age of the order 
of 100 Myr for these most metal-deficient galaxies. The 
chemical evidence for a young age of galaxies with Z ^ 
Zq /20 is supported by photometric and spectroscopic evi- 
dence. Thuan, Izotov & Lipovetsky (1997) and Papaderos 
et al. (1998) have argued, on the basis of colour pro- 
files and spectral synthesis studies, that the second most 
metal-deficient galaxy known after I Zw 18, the BCD SBS 
0335-052 with Z ~ Z Q /A0, is a young galaxy with age 
less than 100 Myr. Using the same techniques in addi- 
tion to colour-magnitude diagram studies, Thuan, Izotov 
& Foltz (1999a) have shown that the age of the BCD SBS 
1415+437 with Z ~ Z©/21 is also less than 100 Myr. 

In view of the contradictory conclusions reached by 
different authors based on the same HST data set, we 
have decided to address anew the issue of the age of I 
Zw 18 by reexamining the HST data. We use the archival 
HST WFPC2 B, V and R images to construct colour- 
magnitude diagrams for I Zw 18. We complement the 
imaging data with high signal - to - noise ratio Multi- 
ple Mirror Telescope (MMT) spectroscopic observations. 
These data are interpreted with the help of spectral en- 
ergy distributions (SEDs). Spectroscopic observations are 
crucial for this type of analysis as they allow to correct 
SEDs for contamination by emission from ionized gas. 
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Fig. 1. Hubble Space Telescope WFPC2 exposure of I 
Zw 18 in the R band. The overlayed contours delineate 
the central Ha shells surrounding the NW and SE star- 
forming knots. 

Table 1. Integrated photometric properties of I Zw 18 



Parameter 


Main body" 


C component" 


Ref. 


RA(J2000) 


09 h 34 m 02.10 s 


09 h 33 m 59.68 s 


1 


DEC(J2000) 


+55°14'23"8 


+55°14'42'.'3 


1 


V (mag) 


16.03 


19.06 


1 


U - B (mag) 


-0.96 




1,2,3 


B — V (mag) 


-0.21 


-0.08 


1 


V - R (mag) 


0.24 


0.15 


1 


V — I (mag) 


-0.12 




1,2,3 



References: 1 - Dufour et al. (1996b); 2 - Hunter & Thronson 
(1995); 3 - this paper. 

"Magnitudes and colours are corrected for interstellar extinc- 
tion (A v = 0.25 mag, E{B - V) = 0.07 mag, E(V - R) = 0.04 
mag, Dufour et al. 1996b). B, V, R magnitudes are taken from 
Dufour et al. (1996b), U, I magnitudes are measured in this 
paper within the circular aperture with radius r = 10" using 
HST observations by Hunter & Thronson (1995). 

In Sect. 2 we discuss the photometric and spectroscopic 
data. A determination of the distance to I Zw 18 and con- 
straints on its age from (B — V) vs. V and (V — R) vs. R 
colour-magnitude diagrams are given in Sect. 3. In Sect. 4 
age constraints inferred from the spectral energy distribu- 
tions are discussed. Hydrodynamical age constraints from 
observations of large expanding shells of ionized gas are 
considered in Sect. 5. We summarize our results in Sect. 
6. 



2. Observational data 

2.1. Photometric data 

Our main goal here is to study the evolutionary status of 
I Zw 18, the integrated photometric properties of which 
are shown in Table [l], as a whole (i.e. including both the 
main body and the C component) by means of colour- 
magnitude diagrams (CMD). To construct the CMDs, we 
use archival deep B (F450W), V (F555W) and R (F702W) 
HST observations ( PI: Dufour, GO-5434, November 1994) 
where both the main body and the C component were 
imaged with the WFPC2 camera. Additionally, for the 
measurement of the total magnitudes (Table Q), we use U 
(F336W) and I (F814W) PCI images of the main body 
(PI: Hunter, GO-5309, November 1994). U and I observa- 
tions of the C component have not yet been done by HST. 
Images reduced by the standard pipeline at the STScI were 
retrieved from the HST archives. They are listed in Table 

I 

The next steps of data reduction included combining 
all observations with the same filter, removal of cosmic 
rays using the IRAF0 routine CRREJ, and correction for 
geometric distortion by producing mosaic images. The lat- 
ter procedure introduces small photometric corrections of 
~ 0.02 mag for objects near the edge of CCD chip (Holtz- 
man et al. 1995a), like the C component, and are even 
smaller for objects located far from the edge, like the main 
body. 

Figure @ shows the HST WFPC2 image of I Zw 18 in 
the R band. A filamentary low surface brightness (LSB) 
emission pattern extending far away from the main body 
of I Zw 18 can be seen. A number of distinct super-shells 
are visible within this LSB envelope of the main body, 
testimony to a large scale perturbation of the surrounding 
gaseous component and suggesting that emission by ion- 
ized gas contributes a significant fraction of the integrated 
luminosity of the system. This is obviously not the case 
for the C component in which weak Ha emission is barely 
seen and is localized in a compact region (Dufour et al. 
1996b). 

To construct the CMDs, we use the DAOPHOT pack- 
age in IRAF. The point-spread function in each broad- 
band filter is obtained with the PSF routine in an inter- 
active mode, by examining and extracting 5-10 stars 
in relatively uncrowded fields. The photometry of point 
sources is done with a 2-pixel radius circular aperture. In- 
strumental magnitudes are then converted to magnitudes 
within an aperture of radius 0'.'5, adopting corrections of 
-0.19 mag in B and V and -0.21 mag in R (Holtzman et 
al. 1995a). They are finally transformed to magnitudes in 

3 IRAF: the Image Reduction and Analysis Facility is dis- 
tributed by the National Optical Astronomy Observatories, 
which is operated by the Association of Universities for Re- 
search in Astronomy, In. (AURA) under cooperative agreement 
with the National Science Foundation (NSF). 
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Table 2. HST/WFPC2 Archival Data of I Zw 18 



Filter 


WFPC2 quadrant 


PI 


Epoch 


Exposure time 
(sec) 


Image root name 


F336W 


PC 


Hunter 


29 Oct 1994 


4200 


u2cg0101t, u2cg0102t, u2cg0103t 


F450W 


WF3 


Dufour 


03 Nov 1994 


4600 


u2f90102t, u2f90103t 


F555W 


WF3 


Dufour 


03 Nov 1994 


4600 


u2f90104t, u2f90105t 


F702W 


WF3 


Dufour 


03 Nov 1994 


5400 


u2f90101t, u2f90201t, u2f90202t 


F814W 


PC 


Hunter 


30 Oct 1994 


6600 


u2cg0301t, u2cg0302t, u2cg0303t 



the standard Johnson-Cousins UBVRI photometric sys- 
tem using the prescriptions of Holtzman et al. (1995b). 

2.2. Spectroscopic data 

Spectroscopic observations of I Zw 18 were carried out 
with the MMT on the nights of 1997 April 29 and 30. 
Signal-to- noise ratios S/N ib 50 were reached in the contin- 
uum of the bright central part. Observations were made in 
the blue channel of the MMT spectrograph using a highly 
optimized Loral 3072 x 1024 CCD detector. A 1'.'5 x 180" 
slit was used along with a 300 groove mm -1 grating in 
first order and an L-38 second-order blocking filter. This 
yields a spatial resolution along the slit of 0'.'3 pixel -1 , a 
scale perpendicular to the slit of 1.9 A pixel -1 , a spectral 
range 3600 - 7500 A, and a spectral resolution of ~ 7 A ( 
FWHM ). For these observations, CCD rows were binned 
by a factor of 2, yielding a final spatial sampling of $.'6 
pixel -1 . The total exposure time was 180 minutes bro- 
ken up in six subexposures, 30 minutes each, to allow for 
a more effective cosmic-ray removal. All exposures were 
taken at small airmasses ( ^ 1.1 - 1.2 ), so no correction 
was made for atmospheric dispersion. The seeing during 
the observations was 0'.'7 FWHM. The slit was oriented 
in the position angle P.A. = -41° to permit observations 
of the NW and SE components and the SE tip of the C 
component simultaneously. The spectrophotometric stan- 
dard star HZ 44 was observed for flux calibration. Spectra 
of He-Ne-Ar comparison lamps were obtained before and 
after each observation to provide wavelength calibration. 

Data reduction of spectral observations was carried 
out at the NOAO headquarters in Tucson using the IRAF 
software package. This included bias subtraction, cosmic- 
ray removal and flat-field correction using exposures of 
a quartz incandescent lamp. After wavelength mapping, 
night-sky background subtraction, and correcting for at- 
mospheric extinction, each frame was calibrated to ab- 
solute fluxes. One-dimensional spectra were extracted by 
summing, without weighting, different numbers of rows 
along the slit depending on the exact region of interest. 
We have extracted spectra of two regions: (1) the bright- 
est part of the main body with size 1'.'5 x 8'.'5 centered on 
the NW component and (2) the southeastern tip of the C 
component within an aperture l'/5 x 3". The observed and 



extinction-corrected emission-line intensities in the main 
body of I Zw 18 are listed in Table § 

The ionic and elemental abundances have been derived 
following the procedure by Izotov et al. (1994, 1997c). The 
extinction coefficient C(H/3) and the absorption equiva- 
lent width EW(&bs) for the hydrogen lines obtained by 
an iterative procedure arc included in Table || together 
with the observed flux F(Hf3) and the equivalent width 
EW(H{3) of the H/3 emission line. The electron temper- 
ature T e (0 III) was determined from the [OIII] A4363/ 
(A4959 + A5007) flux ratio and the electron number den- 
sity iV e (S II) from the [S II] A6717/A6731 flux ratio. The 
ionic and elemental abundances are shown in Table |] to- 
gether with ionization correction factors (ICFs) . They are 
in good agreement with the abundances derived by Skill- 
man k Kennicutt (1993), Izotov & Thuan (1998), Vi'lchez 
& Iglesias-Paramo (1998) and Izotov et al. (1999). 

3. Stellar population ages from colour-magnitude 
diagrams 

The superior spatial resolution of HST WFPC2 images 
combined with the proximity of I Zw 18 permits to re- 
solve individual bright stars and study stellar populations 
in this galaxy by means of colour - magnitude diagrams. 
Such a study has been done by Hunter & Thronson (1995) 
for the main body in the U, V and I bands. Their pho- 
tometry shows a broad main sequence of massive stars 
and blue and red supergiants. The NW component con- 
tains the brightest and reddest, presumably most evolved 
stars, spanning a range 2-5 Myr in age. The stars in the 
SE component are likely to be even younger. 

Dufour et al. (1996b) were able to resolve the C com- 
ponent into stars and found it to be older than the NW 
and SE components. On the basis of (B — V) vs. V colour- 
magnitude diagrams, they concluded that star formation 
started in the main body several tens of Myr ago and 
in the C component several hundred of Myr ago. In this 
section we re-analyze the observations by Dufour et al. 
(1996b) and re-examine their conclusions concerning the 
age of the stellar content of I Zw 18. In order to compare 
observed to theoretical CMDs and derive ages, a precise 
distance to I Zw 18 is needed, which we discuss next. 
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Table 3. Emission line intensities in the main body of I 
Zw 18 



Table 4. Heavy Element Abundances in the main body 
of I Zw 18 



Ion 



F(A)/F(H/3) /(A)/J(H/J) 



3727 [O II] 

3835 H9 

3868 [Ne III] 

3889 He I + H8 

3968 [Ne III] + H7 

4101 m 

4340 H7 

4363 [O III] 

4471 He I 

4686 He II 

4861 H/3 

4959 [O III] 

5007 [O III] 

5876 He I 

6300 [O I] 

6312 [S III] 

6563 Ha 

6678 He I 

6717 [S II] 

6731 [S II] 

7065 He I 

7135 [Ar III] 

C(H/3) dex 
F(H/3) a 
EW(U/3) A 
EW(abs) A 



0.3473±0.0035 
0.0195±0.0021 
0.1519±0.0028 
0.1191±0.0026 
0.1446±0.0026 
0.2072±0.0027 
0.4204±0.0034 
0.0639±0.0022 
0.0263±0.0019 
0.0249±0.0019 
1.0000±0.0053 
0.6927±0.0042 
2.0835±0.0092 
0.0838±0.0017 
0.0090±0.0014 
0.0091±0.0012 
3.0488±0.0126 
0.0291±0.0012 
0.0334±0.0014 
0.0216±0.0014 
0.0251±0.0011 
0.0195±0.0011 



0.3547±0.0038 
0.0971±0.0141 
0.1537±0.0029 
0.1956±0.0052 
0.2170±0.0048 
0.2720±0.0042 
0.4701±0.0043 
0.0628±0.0022 
0.0257±0.0019 
0.0240±0.0019 
1.0000±0.0056 
0.6605±0.0042 
1.9822±0.0092 
0.0768±0.0017 
0.0081±0.0013 
0.0082±0.0011 
2.7450±0.0129 
0.0259±0.0012 
0.0297±0.0013 
0.0192±0.0013 
0.0221±0.0011 
0.0172±0.0011 



0.090±0.005 
3.50±0.02 
68±5 
3.0±0.1 



Mn units of 10 



ergs s cm 



3.1. The distance to I Zw 18 

For the nearest galaxies, distances can be derived from 
colour-magnitude diagrams themselves by measuring the 
apparent magnitude of the tip of the red giant branch 
clump (e.g. Schultc-Ladbeck, Crone & Hopp 1998; Lynds 
et al. 1998). However, red giants are too faint to be seen in 
more distant galaxies (such as I Zw 18), and other meth- 
ods must be used. Any such distance determination should 
always be checked for consistency, i.e. galaxy properties 
derived from CMD analysis such as the age of the stel- 
lar populations or the luminosities of the brightest stars 
should be compatible with other known observed charac- 
teristics of the galaxy. 

A distance of 10 Mpc to I Zw 18 has generally been 
adopted by previous authors (Hunter & Thronson 1995, 
Dufour ct al. 1996b and Aloisi et al. 1999). This assumes 
that the observed heliocentric radial velocity of the galaxy 
~ 740 km s _1 is pure Hubble flow velocity and a Hubble 
constant H = 75 km s _1 Mpc -1 . Adopting this distance 
would lead to a conflict with the well-observed ionization 
state of I Zw 18. At 10 Mpc the brightest stars observed 
in the main body and in the C component have absolute 
V magnitudes fainter than -8 and -6 mag respectively 
(Hunter & Thronson 1995; Dufour ct al. 1996b; Aloisi et 



Parameter 



Value 



T e (0 III)(K) 
T e (0 II) (K) 
T e (S III)(K) 
iV e (S II)(cm- 3 ) 

O+/H+(xl0 5 ) 
O++/H+(xl0 5 ) 
O +3 /H+(xl0 5 ) 
O/H(xl0 5 ) 
12 + log(0/H) 

Ne++/H+(xl0 5 ) 

ICF(Ne) 

log(Ne/0) 

S+/H+(xl0 7 ) 
S++/H+(xl0 7 ) 
ICF(S) 
log(S/0) 

Ar++/H+(xl0 7 ) 

ICF(Ar) 

log(Ar/Q) 



19300±400 
15500±300 
17700±300 
10±10 

0.280±0.014 
1.221±0.058 
0.042±0.004 
1.544±0.060 
7.189±0.017 

0.194±0.010 
1.26 
-0.799±0.033 

0.457±0.022 
2.619±0.377 
1.67 
-1.478±0.045 

0.473±0.030 
2.00 
-2.213±0.026 



al. 1999). If that is the case, comparison with evolution- 
ary tracks implies that the most massive stars in the main 
body and C component would have masses less than 15 
Mq and 9 Mq respectively (Dufour et al. 1996b; Aloisi et 
al. 1999). According to Vacca (1994), O stars have masses 
exceeding 13 Mq. Thus, only very late O stars would be 
present in the main body. This conclusion is in severe con- 
tradiction with the observed ionization state of I Zw 18. 
Indeed, the equivalent width of the H/3 emission line ex- 
pected from a stellar population having an upper mass 
limit of 15 Mq is < 10 A, while the observed H/3 equiva- 
lent widths in the NW and SE components lie in the range 
60 - 130 A, implying the presence of stars with masses 
40 - 50 Mq. If the upper stellar mass limit of 9 Mq de- 
rived by Dufour ct al. (1996b) and Aloisi et al. (1999) for 
the C component is correct, then ionized gas should not 
be present in this component because of the lack of O 
and early B stars. But Ha and H/3 are clearly observed 
(Dufour et al. 1996ab; Izotov & Thuan 1998; van Zee et 
al. 1998; this paper) in the C component. Izotov & Thuan 
(1998) derived EW(H(3) = 6 A, which, after correction for 
underlying stellar absorption, would result in a value as 
high as 10 A. Thus, late O and early B stars with masses 
as high as 15 M© must be postulated in the C component. 
We argue therefore that the stellar absolute magnitudes 
derived by Dufour ct al. (1996b) and Aloisi et al. (1999) 
from their CMDs are too faint because they are based on 
too small an adopted distance. 
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Fig. 2. Colours of point sources in the main body as a 
function of angular distance from the southeastern tip of 
the SE component, for different upper limits of the photo- 
metric uncertainty a. The boundaries of the SE and NW 
components are shown by vertical lines. Faint red stars 
with i? ^2 25 mag and (V — R) iS 0.4 mag are shown by 
open circles. 

The real distance to I Zw 18 appears to be consider- 
ably larger. Correction of the I Zw 18 heliocentric radial 
velocity to the centroid of the Local Group and for a Vir- 
gocentric infall motion of 300 km s _1 ( Kraan-Korteweg 
1986 ) gives a velocity of 1114 km s _1 . If the Hubble con- 
stant is in the currently accepted range of 55 - 70 km s _1 
Mpc" 1 this would give a distance between ~ 16 and ~ 20 
Mpc for I Zw 18. 

Because of uncertainties introduced by the statistical 
nature of the Virgocentric infall correction and those in 
the value of the Hubble constant, we prefer to base our 
estimate of the distance to I Zw 18 on two firm obser- 
vational results: 1) the presence of ionized gas in the C 
component and 2) the presence of WR stars in the NW 
component. Concerning the ionized gas in the C compo- 
nent, the distance should be increased to a value of ~ 20 
Mpc. Increasing the distance by a factor of 2 would make 
the most massive stars more luminous by a factor of 4 



and push the mass upper limit to ~ 15 Mq. These more 
massive stars would then provide enough ionizing photons 
to account for the observed emission lines in the C com- 
ponent. Concerning the WR stars, they have been seen 
both spectroscopically (Izotov et al. 1997a; Legrand et al. 
1997) and photometrically (Hunter & Thronson 1995; de 
Mello et al. 1998) in the NW component of I Zw 18, in 
the region where the brightest post-main-sequence stars 
are located. The existence of WR stars implies the pres- 
ence of very massive stars in the NW component. De Mello 
et al. (1998) using Geneva stellar evolutionary tracks for 
massive stars with enhanced stellar wind and with heavy 
element mass fraction Z — 0.0004 found the minimum 
initial mass M m j n for stars evolving to WR stars to be 
~ 90 Mq and that the WR stage is very short-lived, be- 
ing only ~ 0.8 Myr, if the instantaneous burst model is 
adopted. The models used by de Mello et al. (1998) do not 
take rotation into account, which may decrease M m i ni but 
probably by not more than a factor of 1.5 (Langer & Heger 
1998; Meynet 1998). Thus the observation of WR stars in 
the NW component implies that post-main-sequence mas- 
sive stars with M 40 - 60 Mq and with lifetimes ^ 3 - 
5 Myr (Fagotto et al. 1994; Meynet et al. 1994) must be 
present in I Zw 18. This short time scale is in excellent 
agreement with the age of 5 Myr derived from the equiv- 
alent width of H/3 in the main body (Table g), using the 
calibration by Schaerer & Vacca (1998). To accomodate 
this small time scale, we are forced again to increase the 
distance of I Zw 18 to ~ 20 Mpc. Increasing the distance 
would increase the luminosity of all stars and the location 
of the brightest stars in the CMD of the main body can 
be accounted for by isochrones with age as short as ~ 5 
Myr. 

In summary, three different lines of argument — a non 
Hubble flow velocity component, the ionization state of I 
Zw 18 and the presence of WR stars in its NW component 
— have led us to conclude that I Zw 18 is twice as distant 
as thought before. We shall thus adopt a distance of 20 
Mpc to I Zw 18. At this distance, 1" = 97 pc. 

3.2. The resolved component and colour-magnitude 
diagrams 

In Fig. U we show the spatial distribution of the (V — R) 
colours of point sources along the main body for increas- 
ing photometric uncertainties a. Panels (a) to (d) show 
CMD data points with successively larger photometric er- 
rors, but smaller than the a value given in each panel. The 
origin is taken to be at the south-eastern tip of the main 
body. The boundaries of the SE and NW components are 
shown by vertical lines. In the most crowded fields of the 
SE and NW starburst components, V and R magnitudes 
are measured with a precision better than 0.15 mag only 
for very few stars (Fig. ||a). However, if the upper lim- 
its are greater than > 0.25 mag, then almost all stars in 
these two components are recovered as suggested by the 
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Fig. 3. Photometric uncertainties a in the determination 
(left) and in the main body (right) of I Zw 18. Note the 
as compared to the C component. The larger errors are 
crowding. 

similarity of the distributions of stars in panels (c) and 
(d) . Figure shows the distribution of photometric errors 
as a function of apparent B, V and R stellar magnitudes 
for the C component (Figs. ||a - ^p) and the main body 
(Figs. ||d - ||f). The photometric errors in the V band at 
a fixed magnitude are similar to those derived by Hunter 
& Thronson (1995). In the main body, the errors remain 
rather large even for the brightest stars ( a ~ 0.2 mag at 

V ~ 22 mag ). This is due to the spatially varying back- 
ground in the highly crowded region of the main body. 
These photometric errors are especially large at the faint 
magnitudes of red evolved stars ( a ~ 0.4 mag at V <~ 
28 mag ). By contrast, the photometric errors in all filters 
are lower for the sources in the C component, being ~ 
0.1 mag at V ~ 25 mag and increasing to ~ 0.3 mag at 

V ~ 28 mag, because stellar crowding and contamination 
of the stellar background by the ionized gas emission are 
smaller. 

In the following CMD analysis, we shall consider only 
point sources which are brighter than 27 mag in each band 
and those for which the photometric errors do not exceed 




22 24 26 28 



Magnitude 

of apparent magnitudes of point sources in the C component 
much larger errors at a fixed magnitude in the main body 
due to the contamination of ionized gas and severe stellar 



0.25 mag. For the C component, each of these two condi- 
tions selects out nearly the same point sources, as 27 mag 
sources have a mean error of ~ 0.25 mag. However, in the 
main body, some bright sources are rejected because they 
have photometric errors larger than 0.25 mag (Fig. |||). 

The correction for internal extinction poses a prob- 
lem. The extinction for the main body can be estimated 
from the optical emission-line spectrum. However, differ- 
ent measurements give somewhat different values for the 
extinction. Dufour et al. (1996b) have adopted a value 
Ay — 0.25 mag based on the spectroscopic observations 
of Skillman & Kennicutt (1993). The value of C(H/3) = 
0.09 derived in this paper (Table ||) for the main body 
translates into a mean Ay = 0.19 mag. The extinction 
may differ in the NW and SE components and it is not 
possible to estimate it in the C component because Balmer 
hydrogen emission lines in its spectrum are not detected 
with high enough signal-to-noise ratio. 

It is worth noting that recent studies reveal that dust 
formation may proceed very efficiently in a low-metallicity 
starburst environment and cause a highly inhomogeneous 
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Fig. 4. {B — V) vs. V colour- magnitude diagram for point 
sources in I Zw 18. Theoretical isochrones from Bertelli et 
al. (1994) for a stellar population with a heavy element 
mass fraction Z = 0.0004 are shown by solid lines labeled 
by the logarithm of age in yr and, in parentheses by the 
maximum stellar mass corresponding to each isochrone. A 
distance of 20 Mpc is adopted. 



extinction pattern. Thuan, Sauvage & Madden (1999b) 
have analyzed ISO observations of the second most metal 
- deficient galaxy known, SBS 0335-052, and found that 
an intrinsic extinction of up to Ay = 19-21 mag is re- 
quired around some stellar complexes to account for the 
galaxy's spectral energy distribution in the mid-infrared. 
The mean extinction in SBS 0335-052 derived from the 
optical spectroscopic observations is much lower ( Ay ~ 
0.5 mag, Izotov et al. 1997b), although it varies signifi- 
cantly along the slit. Thuan et al. (1999b) found that as 
much as 3/4 of the massive stars in SBS 0335-052 may be 
hidden by dust. 

Keeping in mind these findings, we adopt for simplic- 
ity and for lack of more information a spatially constant 
extinction of Ay — 0.19 mag, E(B — V) = 0.06 mag and 
E(V -R) = 0.04 mag for the main body of I Zw 18 (Table 
||), while the extinction for the C component is taken to be 
zero. An inhomogeneous and locally much higher extinc- 



Fig. 5. (V — R) vs. R colour-magnitude diagram for point 
sources in I Zw 18. Theoretical isochrones from Bertelli et 
al. (1994) for a stellar population with a heavy element 
mass fraction of Z = 0.0004 are shown by solid curves 
and they labeled as in Fig. ||. A distance of 20 Mpc is 
adopted. 



tion than the value adopted above will evidently result in 
a larger dispersion in the CMDs and in increased colour 
indices, i.e. in an overestimate of the age for I Zw 18. 

We show the (B — V) vs. V diagram for the C compo- 
nent in Fig. [|a and that for the main body in Fig. ||b. Solid 
lines represent theoretical isochrones (Bertelli et al. 1994), 
for a heavy element mass fraction Z = 0.0004, labeled ac- 
cording to the logarithm of age. The largest stellar mass 
associated with each isochrone is also given in parenthe- 
ses. It may be seen that there exists a well-defined main 
sequence in the C component with a turn-off indicative of 
an age of ~ 15 Myr. There is a very bright point source 
with V ~ 22 mag, which may be in fact a compact star 
cluster as suggested by Dufour et al. (1996b). Although 
no red sources with (B — V) iS 1 are seen, there is a group 
of relatively redder points with 0.2 ^ (B — V) ^ 0.6 in 
Fig. |4^l. Because of their faint magnitudes ( V ~ 26 ), 
some of the red sources may be attributed to photometric 
uncertainties ( a ~ 0.2 mag at B and V ~ 26, resulting 
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in cr(B — V) ~ 0.3 mag, Fig. || ). The fact that there are 
also points that scatter to the blue at that faint magni- 
tude would support that hypothesis. However, even if we 
do accept that all these faint red stars are real, compari- 
son with theoretical isochrones says that their age cannot 
exceed 100 Myr. 

As for the main body, the main sequence turn-off at V 
~ 22 mag implies an age of ~ 5 Myr ( Fig. ||b ). Red stars ( 
(B—V) ^ 0.2 ) are seen to possess a large range of absolute 
magnitudes ( ~ 3 — 4 mag ) implying that star formation 
has been undergoing in the main body in the last ~ 15 
- 30 Myr. This conclusion is essentially the same as that 
by Dufour et al. (1996b) except for differences introduced 
by the larger distance to I Zw 18 which results in smaller 
ages. We note that the spread of the points in the (B — V) 
vs. V diagram of the main body (Fig. ||b) is larger than 
the one in the CMD of the C component. It is probably 
not only due to evolutionary effects, but also to larger 
photometric uncertainties in the main body because of 
more crowding and contamination from gaseous emission 
as discussed above (cf. Fig. ||). 

(V — R) vs. R CMD diagrams for the C component 
and the main body along with theoretical isochrones from 
Bertelli et al. (1994) for Z = 0.0004 (solid lines) are shown 
in Fig. ||. As in Fig. ||, the logarithm of age and the max- 
imum stellar mass are given for each isochrone. 

Figure ||a shows a well-defined main-sequence for com- 
ponent C corresponding to an age of ~ 15 Myr. Several red 
and faint (V J; 25 mag) stars with (V — R) Ss 1.1 mag are 
seen. They probably are red supergiants and/or massive 
young luminous asymptotic giant branch (AGB) stars, ~ 
2 mag brighter than the older AGB stars observed in the 
BCD VII Zw 403 (Lynds et al. 1998; Schulte-Ladbeck et 
al. 1998) and Local Group galaxies ( e.g. Gallart, Aparicio 
& Vilchez 1996 ). The location of these red stars in the 
CMD suggests an age of ^ 100 Myr. We conclude that 
star formation in the C component probably began ~ 100 
Myr ago and finished ~ 15 Myr ago. 

The spread of points in the main body ( Fig. ||b ) 
is similar to that in the {B - V) vs. V CMD (Fig. §). 
Again, red stars in the main body span a range by 3 - 4 
magnitudes in brightness suggesting that star formation 
has occurred during the last ~ 15 - 30 Myr. The age of ~ 
5 Myr derived for the brightest stars in the main body is 
in agreement with that inferred from spectral observations 
( see Sect. 4 ). 

Our derived upper age limit of ~ 100 Myr for I Zw 18 
differs by a whole order of magnitude from the age of up to 
1 Gyr derived by Aloisi et al. (1999) in their CMD analysis 
of the archival B,V, I HST data. This large age difference 
comes mostly from the doubling of the distance to I Zw 
18. Aloisi et al. (1999) found an appreciable number of 
faint red sources with V = 26 - 27 mag and (V — I) 
= 1-2 mag. With our new distance of 20 Mpc, these 
sources have absolute I magnitude of -6.5 mag, or 2 mag 
brighter than the absolute / magnitude of older and fainter 



AGB stars detected in the BCD VII Zw 403 (Lynds et al. 
1998; Schulte-Ladbeck et al. 1998). Because of their high 
luminosities, these red stars cannot be interpreted as old 
AGB stars as proposed by Aloisi et al. (1999). Rather, 
it is more likely that they are young stars. At the new 
distance, the location in the (V — I) vs. V diagram of the 
faint red stars seen by Aloisi et al. (1999) is consistent 
with an age of ~ 50 Myr if the Geneva tracks are used, 
and of Si 100 Myr if the Padova tracks are used. These red 
stars with R = 25 - 26 mag and (V — R) ~ 0.6 are also 
seen by us in the (V -R)vsR CMD (Fig. |b) and their 
age is again consistent with £5 50 Myr as inferred from 
the Padova isochrones. Aloisi et al. (1999) restricted their 
CMDs to points with photometric errors smaller than a 
= 0.2 mag. With such a cut, those authors found that 
the faint red stars are located mainly in the southeastern 
part of the main body. Adopting the same cut, we come 
to the same conclusion as can be seen in Figure ^|b, where 
these stars are shown by open circles. They are seen to be 
located mainly in the relatively uncrowded region between 
the SE and NW components. However, if a slightly larger 
cut of a = 0.25 mag is used instead, the distribution of 
faint red stars in the main body becomes more uniform 
(Figure ||c). The ages of the faint red stars of ~ 100 Myr 
derived here, should be considered as upper limits. If these 
sources are subject to local high extinction, their ages will 
be decreased. 

In summary, the upper ages derived for point sources 
from colour-magnitude diagrams in both the main body 
and C component of I Zw 18 do not exceed 100 Myr. 

Figure || shows the radial colour distribution of point 
sources along the body of the C component with the origin 
taken to be at the southeastern tip. While the (B — V) 
colour does not show a gradient (Fig. ^|a), there is a weak 
trend for the (V — R) colour of stars to become redder 
away from the southeastern tip towards the northwestern 
tip of the C component. This trend was also discussed by 
Dufour et al. (1996b) and may reflect propagating star 
formation in the C component from the NW to the SE. 

The dots in Fig. || show the colours of the unresolved 
stellar continuum of the C component down to a surface 
brightness limit of 25 R mag arcsec -2 , averaged within 
regions of size 0'.'4 x 0'.'4. The mean colours of the diffuse 
component were determined to be <B — V> = +0.05 mag 
and < V — R> = +0.16 mag with a standard deviation of 
0.01 mag around the mean. 

A colour gradient in the diffuse stellar continuum is 
more evident in (B — V) than in (V — R), amounting to 
(0.15+0.03) mag kpc^ 1 and (0.05+0.025) mag kpc" 1 , re- 
spectively. The diffuse stellar continuum distribution fol- 
lows very closely that of the resolved point sources, al- 
though a small difference in the (V — R) colour of the 
resolved and unresolved components exists in the south- 
eastern part. These results imply that, contrary to the 
majority of BCDs (e.g. Papaderos et al. 1996a), in which 
star-forming regions are immersed in an extended and 
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Fig. 6. Colours of point sources (filled circles) and of the 
diffuse unresolved low-surface-brightness (LSB) emission 
(dots) in the C component as function of angular distance 
from the southeastern tip of the C component. A slight 
(V — R) colour gradient for resolved sources is present 
which suggests that bright stars in the northwestern part 
of the C component are slightly redder compared to stars 
in the southeastern part. 



much older stellar LSB envelope, the blue unresolved stel- 
lar continuum in the C component is of comparable age 
or even formed coevally with the population of resolved 
sources discussed in this section. The integrated photo- 
metric properties of the unresolved components in both 
the main body and C component are listed in Table ^|. 
They are derived from the total luminosity of each com- 
ponent after the summed up emission of point-like sources 
has been subtracted. Quantities listed in Table [| referring 
to the main body of I Zw 18 are corrected for intrinsic 
extinction. 

Note that while the emission of the resolved compo- 
nent in the main body of I Zw 18 is contributed by stellar 
sources, the unresolved emission includes both stellar and 



Table 5. Photometric properties of resolved and diffuse 
emission in I Zw 18 

Parameter Main body" C component 6 



a) Total 

V (mag) 16.05 19.20 
B — V (mag) -0.18 0.00 

V — R (mag) 0.20 0.14 

b) Resolved population 

V (mag) 18.35 20.69 
B — V (mag) -0.41 -0.03 

V - R (mag) 0.13 0.11 

c) Unresolved population 

V (mag) 16.16 19.46 
B — V (mag) -0.16 0.01 

V - R (mag) 0.21 0.16 



a Corrected for extinction as derived from spectroscopic obser- 
vations (A v = 0.19 mag, E(B - V) = 0.06 mag, E(V - R) = 
0.04 mag). The measurements have been done within a circular 
aperture of 10" in radius. 
b Not corrected for extinction. 

gaseous contributions. It follows from Table || that most 
of the light comes from the unresolved stellar component. 
The contribution of resolved sources to the total light of 
both the main body and C component is 25%. 

4. Synthetic spectral energy distribution 

As discussed before, the stellar emission in the main body 
of I Zw 18 is strongly contaminated by emission of ion- 
ized gas from supergiant HII regions. Therefore, to derive 
the BCD's stellar populations age, synthetic spectral en- 
ergy distributions which include both stellar and ionized 
gaseous emission need to be constructed. By contrast, the 
contamination of the light of the C component by gaseous 
emission is small, and photometric and spectral data give 
us direct information on its stellar populations. 

To analyze the stellar populations in the young ion- 
izing clusters of the main body, we use SEDs calculated 
by Schaerer (1998, private communication) for a heavy 
element mass fraction of Z = Z Q /20 (SEDs with lower 
metallicity are not available) and ages in the range of t = 
2-10 Myr. As for the C component, the absence of strong 
ionized gas emission implies that its stellar population is 
older than 10 Myr. We have therefore calculated a grid 
of SEDs for stellar populations with ages between 4 Myr 
and 20 Gyr and heavy element abundances Z/Z Q = 10~ 2 , 
using the stellar isochrones of Bertelli et al. (1994) and 
the compilation of stellar atmosphere models of Lejeune, 
Cuisinier & Buser (1998). An initial mass function (IMF) 
with a Salpeter slope equal to -2.35, an upper mass limit 
of 100 Mq and a lower mass limit of 0.6 M Q are adopted. 
The observed gaseous spectral energy distribution is then 
added to the calculated stellar spectral energy distribu- 
tion, its contribution being determined by the ratio of the 
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observed equivalent width of the H/3 emission line to the 
one expected for pure gaseous emission. To calculate the 
gaseous continuum spectral energy distribution, the ob- 
served H/3 flux and the electron temperature have been 
derived from the spectrum at each point along the slit. 
The contribution of bound-free, free-free and two-photon 
continuum emission has been taken into account for the 
spectral range from to 5 /xm (Aller 1984; Ferland 1980). 
Emission lines are superposed on top of the gaseous con- 
tinuum SED with intensities derived from spectra in the 
spectral range A3700 - 7500 A. Outside this range, the 
intensities of emission lines (mainly hydrogen lines) have 
been calculated from the extinction-corrected intensity of 
H/3. 

4--1- The main body 

Fig. |?j shows the observed spectrum of the main body 
( NW + SE components ) along with the synthesized 
gaseous and stellar continuum for a composite stellar pop- 
ulation of ages 2 and 5 Myr contributing respectively 40 
% and 60 % of the total mass, in the approximation of an 
instantaneous burst of star formation. The line intensities, 
electron temperature, electron number density and heavy 
element abundances are taken from Tables || and |]. Com- 
parison with synthesized SEDs of different ages shows that 
there is good agreement with the observed SED only for 
an age of 5 Myr. However, the synthesized SED is system- 
atically below the observed SED in the blue part of the 
spectrum with A ^ 4000 A. Hence, an additional younger 
stellar population with age ~ 2 Myr is required. We con- 
clude that the emission of the main body is dominated by 
a very young stellar population with age of ^ 5 Myr. This 
age is considerably shorter than the one of 15 - 20 Myr 
derived by Aloisi et al. (1999) for their postulated second 
burst of star formation, and much shorter still than the 
30 Myr - 1 Gyr time scale for the first episode of star 
formation. 

The contribution of ionized gas to the SED of star- 
forming regions in the main body is not as large in I Zw 
18 as in SBS 0335-052 (Papaderos et al. 1998), as can 
be seen by the equivalent width of H/3 and the observed 
intensities of other emission lines relative to H/3. In the 
brightest part of the main body, the equivalent width of 
the H/3 emission line is only 68 A (Table |]) as compared 
to £ 200 A in SBS 0335-052. 

4-2. The C component 

Fig. H shows the spectrum of the southeastern part of the 
C component. It is fitted very well by a single stellar pop- 
ulation with age ~ 15 Myr. This is in excellent agreement 
with the age derived earlier from the main-sequence turn- 
off in the CMD of the C component, supporting our large 
adopted distance of 20 Mpc to I Zw 18. To illustrate the 
sensitivity of synthesized SEDs to age determination, we 
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Fig. 7. MMT spectrum of the main body of I Zw 18 on 
which is superposed a synthetic continuum including both 
gaseous and stellar emission for a composite stellar popu- 
lation composed of a young population with age 2 Myr (40 
% of total mass) and an older one with age 5 Myr (60 % 
of total mass) (thick line). The observed spectral energy 
distribution is corrected for an extinction Ay = 0.19 mag, 
derived from the Balmer hydrogen emission lines decre- 
ment. 

have also shown in Fig. || the SED of a stellar popula- 
tion with age 40 Myr, the value adopted by Dufour et 
al. (1996b). The agreement is not so good, the synthe- 
sized fluxes being systematically smaller than the observed 
fluxes for A < 5000 A. 

We conclude that evolutionary synthesis models fur- 
ther constrain the age of the C component to the range 
15 - 40 Myr. 

4-3. A 10 Gyr old stellar population? 

Although there is compelling evidence that the stellar 
emission in I Zw 18 ( in both the main body and the C 
component ) is due to stellar populations with age 5s 100 
Myr, the presence of a very faint older stellar population 
with age ~ 10 Gyr cannot be definitely ruled out, since 
such a population would not be visible in CMDs with V 
5a 27.5 mag and would cause a non-detectable effect in 
photometric and spectroscopic data. 

Figure |] shows composite SEDs resulting from a mix- 
ture of a young 15 Myr and an old 10 Gyr stellar popu- 
lation. Each SED is labeled by the mass fraction of the 
old stellar population. The SED with comparable masses 
of young and old stellar populations ( labeled as 0.5 ) is 
indistinguishable from that of a pure young stellar popula- 
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Fig. 8. Spectrum of the C component of I Zw 18 (thin 
line) on which are superposed model stellar population 
SEDs with ages 15 Myr and 40 Myr (thick lines). The 
spectrum is uncorrected for extinction. 

tion ( labeled ). However, SEDs of composite populations 
where the mass of the old stellar population is 10 times 
greater than the mass of the younger stellar component 
can be excluded as they are systematically redder than 
the observed spectrum for A ^ 5500 A. 

We cannot exclude therefore an underlying very old 
stellar population with mass comparable to that of the 
young stellar population in the C component. It must be 
however spatially coincident with the young stellar pop- 
ulation because the colours are nearly constant. This is 
very unlikely as photometric studies of other BCDs (e.g. 
Papaderos et al. 1996ab) show redder colours towards the 
outer parts implying that the old stellar population is spa- 
tially distinct and more extended than the young stellar 
population. 

In summary, there is no need for a stellar popula- 
tion older than ~ 100 Myr to account for photomet- 
ric and spectrophotometric properties of I Zw 18. Deep 
near - infrared photometric observations are needed to put 
stronger constraints on a older stellar population. 

5. Age constraints from the ionized gas shell 
structures 

Many shell structures are seen in Ha images of I Zw 
18 (Fig. Q), produced by the combined effect of stellar 
winds from massive stars and supernovae. Knowledge of 
the radii, the number of massive stars, the ambient den- 
sity, the geometry of the gas distribution and the prop- 
erties of stellar winds from low-metallicity massive stars 



Fig. 9. Spectrum of the C component of I Zw 18 on which 
are superposed a modeled SED for a composite stellar 
population including a young stellar population with age 
15 Myr and an older one with age 10 Gyr. Each SED 
is labeled by the relative mass fraction of the old stellar 
population. No interstellar extinction correction has been 
applied. 

and of supernovae allows, in principle, to estimate the age 
of the stellar populations responsible for these large-scale 
structures. However, this knowledge is still rather uncer- 
tain and stellar population age estimates from shell struc- 
tures should be considered to be more qualitative than 
quantitative. 

If the largest shells were produced by the evolution 
of the NW and SE components, then the shell structures 
would have radii <J 1200 pc. However, the centers of sym- 
metry for the large shells are not coincident with the NW 
component. It is more likely that the large shells are pro- 
duced by older stellar clusters as suggested by Dufour et 
al. (1996b), and in that case, the structures would have 
radii of - 400 - 800 pc. 

The radius of a superbubble produced by a population 
of stars with stellar winds and by supernovae is given, 
respectively, by (McCray & Kafatos 1987): 

R s = 269 pc {N w L 38 /n Q )^ 5 t 3 7 /5 , (1) 
R s = 97 pc (NsnE^/ho) 1 / 5 ^ 5 , (2) 

where Rs is in pc, Nw is the number of stars with stellar 
wind, L 38 = Lw/{W 38 ergs s -1 ), Lw being the mechani- 
cal luminosity of a single stellar wind, tr = i/(10 7 yr), no 
is the density of the ambient gas, N$n is the number of 
supernovae, E 5 i — Esn/(10 51 ergs), Esn being the me- 
chanical energy of a single supernovae. 
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Consider first the stellar wind mechanism. The me- 
chanical energy of the stellar wind of a single star in our 
Galaxy is L 38 ~ 1 (Abbott, Bieging & Churchwell 1981). 
Taking into account the decrease of wind efficiency with 
metallicity as Z - 5 (Maeder & Meynet 1994), we derive 
L38 = 0.14 for I Zw 18's metallicity ( note that Hunter 
& Thronson (1995) used three times as large a value ). 
Guseva, Izotov & Thuan (1999) derived a number of 66 
WR stars in I Zw 18. This number should be considered 
as a lower value of stars with stellar winds, since massive 
O stars can also contribute. The number density of the 
ambient neutral gas no can be estimated from H I column 
densities (van Zee et al. 1998). It varies from ~ 1 cm~ 3 in 
the central part of the main body where the inner shell is 
observed, to ~ 0.05 cm -3 at the distance of ~ 1 kpc from 
the main body where the largest shells are seen. We adopt 
two limiting values no = 0.1 and 1 cm -3 , within the range 
of number densities discussed by Martin (1996). The WR 
stage in an instantaneous burst at Z ~ Z /5O is short, 
t 7 ~ 0.1 (de Mello et al. 1999). Then, the radius of the 
shell produced by stellar winds of WR stars is Rs ~ 110 
pc, if no = 1 cm~ 3 and ~ 170 pc, if no = 0.1 cm~ 3 , in 
good agreement with the observed radius ~ 100 pc of the 
inner shell. However, stellar winds alone cannot explain 
the presence of larger structures. 

Consider next the combined effect of supernovae. In 
the NW component, the formation of structures with radii 
as large as 1200 pc should be accounted for. For simplicity, 
we assume the number of supernovae to be equal to the 
number of O stars. The latter can be estimated from the 
total observed flux of the H/3 emission line. The H/3 flux 
in Table || cannot be used as it has not been corrected for 
aperture effects. Instead, we adopt the aperture-corrected 
H/3 flux from Guseva et al. (1999), giving a number of O 
stars equal to 4800. Adopting — 0.5 for the age of the 
main body, Eq. (||) gives Rs = 550 pc and 350 pc for n 
= 0.1 cm~ 3 and 1 cm -3 respectively. 

To explain the presence of structures with radii as large 
as 1200 pc, the age for the NW component should be as 
high as 18 - 39 Myr, larger than the age of ~ 5 Myr in- 
ferred from ionization constraints and the presence of WR 
stars. We conclude therefore that the stellar clusters in the 
NW and SE components ionizing the gas in the main body 
cannot be responsible for the largest structures, and these 
are due to the action of older clusters and stellar associa- 
tions. The largest shell to the west of the NW component 
with radius 800 pc is likely to be connected with the stellar 
association located at 4" to the NW of the NW compo- 
nent, near its center of symmetry. If we assume that a 
starburst occurred 10 Myr ago, that there were as many 
O stars as in the NW component and no = 0.1 cm' 3 , then 
Rs — 840 pc. Because of the weak dependence of Rs on 
£51 and tf ( Eq. || ) even a burst 10 times weaker with 
an age of 20 Myr can account for the radius of the largest 
shell. Hence, star formation in the main body of I Zw 18 is 
likely to have started ~ 20 Myr ago, propagating generally 



from the NW to the SE direction, and continuing the star 
formation started earlier in the C component. This age 
estimate is in agreement with the one derived by Martin 
(1996) with the use of a more complex model of the gas 
distribution, and of the smaller distance of 10 Mpc. 

6. Conclusions 

Our main goal here is to use observed properties of the 
blue compact dwarf galaxy I Zw 18 to put constraints on 
its age: the high ionization state of the gas and the pres- 
ence of WR stars in the main body, the existence of ion- 
ized gas in the C component, and the colour-magnitude 
diagrams from HST images. We were motivated by the 
study of Izotov & Thuan (1999) who have analyzed the 
C/O and N/O abundance ratios of a sample of the most 
metal-deficient BCDs known, including I Zw 18. Those au- 
thors found that these ratios are constant for galaxies with 
Z ^ Zq/20 with a very small dispersion around the mean. 
This strongly suggests that intermediate-mass stars ( M 
Si 8 Mq ) have not had time to release their carbon and 
primary nitrogen production, establishing an age upper 
limit of ~ 100 Myr for very metal-deficient BCDs. 

The conclusion that galaxies with Z ^ Zq / 20 are 
younger than 100 Myr has been supported by photomet- 
ric and spectroscopic studies of two very metal-deficient 
BCDs, SBS 0335-052 ( Z Q /40, Thuan et al. 1997, Pa- 
paderos et al. 1998) and SBS 1415+437 ( Thuan et al. 
1999a). Here we examine the age evidence for I Zw 18. To 
put constraints on the age of I Zw 18, we have followed 
3 independent lines of investigation: colour-magnitude di- 
agram, spectral synthesis and hydrodynamical age con- 
straints. We have arrived at the following main conclu- 
sions: 

1. The distance to I Zw 18 must be increased by a 
factor of 2 from the previously adopted value of 10 Mpc 
to 20 Mpc. Such a distance is required to have stars bright 
and massive enough in I Zw 18 to account for its high state 
of ionization and the presence of Wolf-Rayet stars in its 
NW component. 

2. (B - V) vs. V and (V - R) vs. R CMD studies with 
the new distance of 20 Mpc give ages derived from the 
main sequence turn-off of I Zw 18 of ~ 15 Myr and ~ 
5 Myr for the C component and main body respectively. 
The location of the resolved luminous red stars with Mr 
~ -6 mag and (V — R) ~ 0.6 -1.0 mag is consistent with 
an age ^ 100 Myr for the C component. The star for- 
mation in this component is likely to have stopped 15 - 
20 Myr ago. As for the main body, CMD analysis implies 
that star formation started 20 - 50 Myr ago and still con- 
tinues nowadays. Analysis of shell structures seen in Ha 
images also suggests that star formation in the main body 
began ~ 20 Myr ago in different locations at the NW side 
and has been propagating mainly in the SE direction. The 
age upper limit of ~ 50 Myr derived for the main body 
is a whole order of magnitude smaller that the one de- 
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rived by Aloisi et al. (1999) from CMD analysis of similar 
HST data. The difference is mainly due to the increase in 
distance of I Zw 18 by a factor of 2. 

3. Fits to the spectral energy distributions give ages of 
~ 5 Myr for the main body and ~ 15-40 Myr for the C 
component. 

In summary, all three lines of investigation ( CMDs, 
the distribution of Ha shells and spectroscopy ) lead to 
the same conclusion, that I Zw 18 did not start to form 
stars until 5s 100 Myr ago. This supports the contention 
of Izotov & Thuan (1999) that all very metal-deficient 
galaxies ( Z ^ Z Q /20 ) are young. 
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